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Abstract We consider the optimal dividend problem for the insurance risk process in 
a general Levy process setting. The objective is to find a strategy which maximizes the 
expected total discounted dividends until the time of ruin. We give sufficient conditions 
under which the optimal strategy is of barrier type. In particular, we show that if the 
Levy density is a completely monotone function, then the optimal dividend strategy is a 
barrier strategy. This approach was inspired by the work of Avram et al. (2007) [Annals 
of Applied Probability 17, 156-180], Loeffen (2008) [Annals of Applied Probability 18, 
1669-1680] and Kyprianou et al. (2010) [Journal of Theoretical Probability 23, 547-564] 
in which the same problem was considered under the spectrally negative Levy processes 
setting. 
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1 The model and problem setting 



Let X = {X t : t > 0} be a real-valued Levy process define on a filtered probability 
space (ft, J 7 , F, V) where F = (J-)t>o is generated by the process X and satisfies the usual 
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conditions of right-continuity and completeness. Denote by P x the probability law of 
X when it starts at x. For notational convenience, we write P = P . Let E x be the 
expectation operator associated with P x with E = E . For 9 e R, k{9) denotes the 
characteristic exponent of X given by 

1 1 f°° 

k(6) = - log E{e l6Xt ) = ia9 - -a 2 9 2 + / [e Sx - 1 - i9xl {{xl<1} )U(dx), 

t ^ J-oo 

where a, a are real constants, and II is a positive measure supported on (— oo, oo) \ {0} 
which satisfies the integrability condition 



/oo 
min{l, x 2 }H(dx) < oo. 
-oo 



The characteristics (a, a 2 , II) are called the Levy triplet of the process and completely 
determines its law; II is called the Levy measure; and a is the Gaussian component of X. 
If U(dx) = n(x)dx, then we call n the Levy density. Such a Levy process is of bounded 
variation if and only if a — and f_ min{|x|, l}U(dx) < oo. It is well known that Levy 
process X is a space-homogeneous strong Markov process. 

For He9 = 0, we can define the Laplace exponent of the process X by 

*(6) = -\ogE(e eXt ) = a6 + -a 2 9 2 + / (e 9x - 1 - 9xl { \ xl<1} )U(dx). (1.1) 

That is, 

E(e eXt ) = e**V\ Re9 = 0, t > 0. 
Using Ito's formula, we find that the infinitesimal generator for X is given by 

1 f°° 
Tg(x) = -a 2 g"(x) + ag'(x) + J [g(x + y) - g(x) - 5 '(x)yl {M<1} ]n(^), (1.2) 

for g <E C 2 with compact support, where 1a is the indicator function of set A. 

When I1{(0, oo)} = 0, i.e., the Levy process X with no positive jumps, is called 
the spectrally negative Levy process. In this case, we recall from Bertoin (1998) and 
Kyprianou (2006) that for each q > 0, there exists a continuous and increasing function 
: R — > [0, oo), called the g-scale function defined in such a way that W®(x) = for 
all x < and on [0, oo) its Laplace transform is given by 



Jo 



e -** W to\x)dx = ^— , 9 > 0(g), 

V(0J - q 



where 0(g) = sup{# > : ^(9) = q} is the right inverse of Smoothness of the scale 
function is related to the smoothness of the underlying paths of the associated process. 
The following facts are taken from Kyprianou et al. (2010) and Chan et al. (2010). It is 
known that if X has paths of bounded variation then, for all q > 0, W^|(o,oo) e C 1 (0 ) °°) 
if and only if II has no atoms. In the case that X has paths of unbounded variation, it is 
known that, for all q > 0, IU (9) |(o,oo) £ ^(0, oo). Moreover, if a > 0, then C^O, oo) may 
be replaced by C 2 (0, oo); and if the Levy measure has a density, then the scale function 
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is always differentiable. In particular, if tc is completely monotone, then W^\( 0oo ) E 
C°°(0,oo). 



For general references on Levy processes, we refer the reader, among others, to Bertoin 
(1998), Sato (1999) and Kyprianou (2006). 

In this paper, we only consider the case that II with density n is absolutely continuous 
with respect to Lebesgue measure. For the Levy process X, we consider the following de 
Finetti's dividend problem. Let £ = {L,\ : t > 0} be a dividend strategy consisting of a 
left-continuous non- negative non- decreasing process adapted to the filtration {J-t}t>a of 
X. Specifically, L\ represents the cumulative dividends paid out up to time t under the 
control £ for an insurance company whose risk process is modelled by X. We define the 
controlled risk process = {Iff : t > 0} by C/| = X t - h\. Let r* = inf{t > : U\ < 0} 
be the ruin time when the dividend payments are taken into account. Define the value 
function of a dividend strategy £ by 



where S > is the discounted rate. The integral is understood pathwise in a Lebesgue- 
Stieltjes sense. 

A dividend strategy is called admissible if L,\ + — L\ < C/f for t < r^. In words, the 
lump sum dividend payment is smaller than the size of the available capitals. Let S be 
the set of all admissible dividend policies. The de Finetti dividend problem consists of 
solving the following stochastic control problem: 



and, if it exists, we want to find a strategy G S such that V^(x) = V*(x) for all x > 0. 

This optimization problem goes back to de Finetti (1957), who considered a discrete 
time random walk with step sizes ±1 and showed that a certain barrier strategy maximizes 
expected discounted dividend payments. Optimal dividend problem has recently gained 
great attention in the actuarial literature. For the diffusion risk process, the optimal prob- 
lem has been studied by many authors including Asmussen et al. (2000), Paulsen (2003), 
and Decamps and Villeneuve (2007). It is well known that under some reasonable as- 
sumptions, the optimality in the diffusion process setting is achieved by a barrier strategy 
(see, for example, Shreve et al. (1984)). The general problem for the Cramer- Lundb erg 
risk model was first solved by Gerber in (1969) via a limit of an associated discrete prob- 
lem. Recently, Azcue and Muler (2005) used the technique of stochastic control theory 
and Hamilton- Jacobi-Bellman (HJB) equation to solve the problem. They also included 
a general reinsurance strategy as another possible control. For the Cramer-Lundberg risk 
model with interest, Yuen et al. (2007) investigated some ruin problems in the presence 
of a constant dividend barrier; and Albrecher and Thonhauser (2008) derived the optimal 
dividend strategy which is again of band type and for exponential claim sizes collapses to 
a barrier strategy. In fact, for various risk models, many results in the literature indicate 
that a band strategy turns out to be optimal among all admissible strategies. For a more 
general risk process, namely the spectrally negative Levy process, Avram et al. (2007) 




V m (x) = sup V^(x), 
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gave a sufficient condition involving the generator of the Levy process for the optimal- 
ly of the barrier strategy; Loeffen (2008) connected the shape of the scale function to 
the existence of an optimal barrier strategy and showed that the optimal strategy is a 
barrier strategy if the Levy measure has a completely monotone density; and Kyprianou 
et al. (2010) further investigated the optimal dividend control problem and showed that 
the problem is solved by a barrier strategy whenever the Levy measure of a spectrally 
negative Levy process has a density which is log-convex. 

Motivated by the work of Avram et al. (2007), Loeffen (2008) and Kyprianou et al. 
(2010) for a spectrally negative Levy process, our objective is to consider the optimal 
dividend problem for a general Levy process (not necessarily spectrally negative). In 
contrast to the previous works, our approach does not rely on the theory of scale function 
of a spectrally negative Levy process. Instead, the derivation of our results requires an 
introduction of a new function and the use of the Wiener-Hopf factorization theory. The 
rest of this paper is organized as follows. In Section 2, we present some preliminary results 
which are derived in Rogers (1983), Bertion and Doney (1994) and Kyprianou (2006). In 
Section 3, we discuss a barrier strategy for dividend payments for the risk model of study. 
In Section 4, we give the main results and their proofs. Finally, Section 5 presents some 
examples and ends with a remark. 



2 Two lemmas on the probability of ruin 



In this section, we present two results on the probability of ruin which will be used later. 

The first result is about the complete monotonicity of the probability of ruin. By 
definition, an infinitely differentiable function / e (0, oo) — > [0, oo) is called completely 
monotone if (— l) n f^ n \x) > for all n = 0, 1,2, Denote by r(q) an exponential 
random variable with mean 1/q which is independent of the process X. For q = 0, r(0) 
is understood to be infinite. Furthermore, let 

X t() = inf X t and X r(q) = sup X t , 

W 0<t<r(q) 0<t<r(q) 

be the infimum and the supremum of the the Levy process X killed at the random time 
r(q), respectively. In the case with q = 0, we always assume that 

/oo 
xl { \ x \>i } n(x)dx > 0. (2.1) 
-oo 



We call the characteristic functions E (exp(aX_ T ^)) and E(exp(aX T ^)) the right and 
left Wiener-Hopf factors of X, respectively. According to Rogers (1983), the Wiener-Hopf 
factors are called mixtures of exponentials if there are probability measures H + and H_ 
on (0, oo] such that 

, xx f°° XHJdX) . v f°° XHJdX) 
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Note that if E exp(aX r ^) is a mixture of exponentials, then the probability distribution 

function of 2Lr(q) nas f° rm 

P_(x) = ! e Xx H_(d\), x<0. 

i(0,oo] 

Similarly, the probability distribution function of X T ^ has the form 

P+{x) = [ (1 - e- Xx )H + (d\), x>0. 

J{0,oo] 

For x > 0, let r = inf{t > : x + X(t) < 0}. Define its Laplace transform and the 
probability of ruin as 

$ q {x) = E(e' qr ), 1>(x) = P(r < oo), 

respectively Then, 

1> 9 (x) = P(X T(q) < -x), 

and 

%l)(x) = lim^ ff (a;) = PiK^ < -x). 

q— >0 

Definition 2.1. (Rogers (1983)) The Levy process X has completely monotone Levy 
density if there exist measures yU+,/i_ on (0, oo) such that 

ir{x) = 1 {X>0} [ e- tx ft + (dt) + 1 {X<0} [ e t >_(dt), (2.2) 

</(0,oo) J(0,oo) 

where 

J t ( 1 | t ) 2 (/ i + + / i ~)( rft ) < °°- 

Theorem 2 of Rogers (1983) states that the jump measure U has a completely mono- 
tone density if and only if the Wiener-Hopf factors of X are mixtures of exponential 
distributions. So, we have the following lemma. 

Lemma 2.1. If the jump measure IT has a completely monotone density, then the func- 
tions i/j q (x) and if)(x) are completely monotone in (0, oo). 

The second result is on Cramer's estimate for ruin probability. Roughly speaking, 
under suitable conditions, the probability of ruin decays exponentially when the initial 
capital becomes larger. The following lemma is extracted from Kyprianou (2006, Theorem 
7.6); see also Bertion and Doney (1994, Theorem). 

Lemma 2.2. Assume that X is a Levy process which does not have monotone paths, for 
which 

(i) (HOP holds; 
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(ii) there exists a R > such that ty(—R) = where is given by U.l\) ; 
(Hi) the support ofU is not lattice ifH(R) < oo. 
Then, 



lim e Rx P(r < oo) = k(0, 0) I R 



d(3 



-i 



--R, 



where the limit is interpreted to be zero if the derivative on the right-hand side is infinite. 
For more details, see Kyprianou (2006, Theorem 7.6). 



3 Barrier strategy 



In this section, we consider a simple barrier strategy for dividend payments. Under a 
barrier strategy, if the controlled surplus reaches the level b, then the overflow will be 
paid as dividends; and if the surplus is less than 6, then no dividends are paid out. Let 
= {L\ : t > 0} be a barrier strategy and U b = {U b : t > 0} be the corresponding 
controlled risk process. Note that U\ = X t — L\ and G 5. Moreover, if Uq G [0,b], 
then the strategy £& corresponds to a reflection of the process X — b at its supremum; 
and if t < u = inf{t > : U\ < 0}, the process L\ can be defined by L\ = and 
L\ = sup s<t [X s — b] V 0. Note that L\ is increasing, continuous and adapted such that 
the support of the measure dL b t is contained in the closure of the set {t : U\ = b}. If 
Uq = x > b, L\ has a jump at t — of size x — b to bring U b back to the level b and a 
similar structure afterward. 



Let Vb(x) denote the dividend- value function if a barrier strategy with level b is applied. 
Then, 

V b (x) = E (j o e ~ 5tdL t\ u o = ■ (3- 1 ) 

The following result shows that Vb(x) as a function of x satisfies the following integro- 
differential equations with certain boundary conditions. Note that Paulsen and Gjessing 
(1997) established a similar result for a very general jump- diffusion process. 

Theorem 3.1. Assume that the process X have Laplace exponent ( GUP and the infinites- 
imal generator Y is given by hi. 2) . Let Vb{x) be bounded and twice continuously differen- 
tiate on (0, b) with a bounded first derivative and with the understanding that we mean 
the right-hand derivatives at x = 0. 

(i) IfVb(x) solves 

TV b (x) = SV b (x), < x < b, 
together with the boundary conditions 

V h (x) = 0, x < 0, 

V b (0) = 0, if a 2 > 0, 

K(b) = 1, 

V b (x) = V b (b) + x - b, x > b, 
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then V b (x) is given by A3. 



Furthermore, let ip b {x) be bounded and twice continuously differentiate on (0,6) with 
a bounded first derivative and with the understanding that we mean the right-hand deriva- 
tives at x = 0. 

(ii) If ip b {x) solves 

Tijj b (x) = 0, < x < b, 
together with the boundary conditions 



ip b (x) = 1, x < 0, 
V>&(0) = 1, if a 2 > 
$(&) = o, 



then ip b (x) = Pxijb < oo). 



Proof. The proof of (i) can be done by using the arguments used in Paulsen and Gjessing 
(1997). If cr 2 > 0, the process starting from immediately has a negative value. Hence, 

V b {0) = 0. Applying Ito's formula to e^'^V^l^) gives 



tAr" 



e -^) Vb{u ^ b) = Vb{x)+ e- 5s (T-5)V b (U b s )ds 

Jo 

l>t/\T b 

~ / e- s 'VM)dL\ + M u 
Jo 

where M t is a martingale. Since Vb(U \) = and the support of the measure dL\ is 
contained in the closure of the set {t : Jj\ = b}, taking expectation on both sides of the 
equality above yields 

E (e- s ^V b (U b tATb )) = V b {x) - E ^ T e- Ss dL b )j . (3.2) 

Therefore, (i) follows by letting t — > oo in (13. 2p . The proof of (ii) is entirely analogous to 
the proof of (i). □ 

The integro-differential equation 

Th(x) = 8h(x), x > 0, 

has, apart from a constant factor, a unique nonnegative solution h(x). This together with 
Theorem 3.1 gives 

f MeI o < x < b 

+™ lb (3 ' 3) 

I x - a + iwy x>k 

In particular, if II{(0, oo)} = 0, then h becomes the £-scale function and (I3.3P reduces 
to Proposition 1 of Avram et al. (2007), which was proved by the excursion theory. See 
also Renaud and Zhou (2007) and Zhou (2005) for an alternative approach. 
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4 Main results and proofs 



Define a barrier level by 

b* = sup{6 > : h\b) < h'{x) for all x > 0}, 
where h'(0) is understood to be the right-hand derivative at 0, and is not necessarily finite. 

We now present the main results of the paper which concerns the optimal barrier 
strategy for a general Levy processes. This is a continuation of the work of Avram 
et al. (2007), Loeffen (2008), and Kyprianou et al. (2010) in which only the spectrally 
negative Levy process was considered. 

Theorem 4.1. Suppose that h belongs to C 1 (0, oo) if a = and J™ mm{\x\, l}H(dx) < 

oo and otherwise belongs to C 2 (0, oo). 7/lim^^oo h'(x) = oo and h(x) is convex on [b*, oo). 
Then, the barrier strategy at b* is an optimal strategy. 

For simplicity, we write the Levy density n as 

f m(x), x>0, 
n(x) = < 

[ 7T 2 (-X), X < 0, 

where tti, n 2 are Levy measures concentrated on (0, oo). 

Theorem 4.2. Ifni and n 2 are completely monotone on (0, oo), then the barrier strategy 
at b* is an optimal strategy. 

Before proving the main results, we present several lemmas which are similar to those 
for spectrally negative Levy process. For 5 > 0, we consider the following second order 
integro-differential equation: 

1 f°° 

-a 2 h"(x) + ah'(x) + / [h(x + y)- h(x) - h'(x)yl {lyl<1} ]U(dy) = 5h(x), x > 0. (4.1) 

2 J-OO 

Set 

p{8) = sup{6 : #(0) = 5}. 

We assume that p(S) > 0. For such a p(5), we denote by P p ^ the exponential tilting of 
the measure P with Radom-Nikodym derivative 

dppiS) 



dP 



= e p(5)X(t)-5t _ 



Under the measure P^ 5 ), the process X is still a Levy process with Laplace exponent 
^p(8) given by 

Let the process X has the Levy triplet (a,a 2 ,n) where a 2 = a 2 , Tl(dx) = n(x)dx = 
e p(5)x ir(x)dx, and 



a = a 



/oo 
( e P(% _ l)yl {mi}7T (y) d y. 
-oo 



Moreover, 

/ e p{5)x ir(x)dx < oo. (4.2) 

J\x\>l 

We refer the reader to Kyprianou (2006) for related discussions. 

Note that the law of X is X under the new probability measure. Let ip(x) be the ruin 
probability for the Levy process X. Then, we have the following result. 

Lemma 4.1. The solutions of equation ( |^.i| ) are proportional to the function (l—ip(x))e p ^ x . 

Proof. We first claim that E(Xi) = ^'(0) is always positive where \I/ is the Levy 
exponent of X. In fact, 

/ — 1 poo 
XTrdx + / x'kdx 
-oo Jl 

/oo /*1 
ye p ^H{y)dy- j yir(y)dy 
-oo J — 1 

since \l/ is strictly convex and increasing on [0,p(5)]. If tt is completely monotone, then it 
follows from Definition 2.1 that it has the form f)2.2p . By using the integrability condition 

we get J p(<5) /i+(^) = 0. Therefore, 
n(x) = e^vr(x) = 1 {X>0} [ e^-^^dt) + l {x<0} [ e^ t+p ^ x ^(dt), 

J(p(S),oo) J{0,oo) 

which shows that tt is completely monotone. So, we see from Lemma 2.1 that if) is also 
completely monotone, and hence ip £ C°°(0, oo). Furthermore, it follows from Theorem 
3.1 (ii) that ip solves 

f^(x) = 0, x > 0, 
together with the boundary conditions 

tf(x) = 1, x < 0, and $(0) = 1 if a 2 > 0, 

where T is the infinitesimal generator for X given by 

1 f°° 
r 9i x ) = lf 2 9'\ x ) + ag'{x) + J [g(x + y) - g(x) - ^(a^l^i}]^)^. 

Then, it can be shown by straightforward calculations that the function 

h(x) = (1 - 4){x))e p{5)x 
satisfies integro-differential equation (14. ip . Hence, the result follows. □ 

Moreover, using arguments similar to those in Avram et al. (2007) and Loeffen (2008), 
we have the following two results. 
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Lemma 4.2. Suppose that h belongs to C 1 (0, oo) if a = and f™ min{|x|, l}H(dx) < oo, 
and otherwise belongs to C 2 (0, oo). If h is also convex on [b*,oo), then 

{T-a)V^(x) <0, for x>b*. 

Lemma 4.3. (Verification lemma) Suppose that £ is an admissible dividend strategy such 
that is twice continuously differentiable and for all x > 

max{(T - a)V ( (x), 1 - V^x)} < 0. 
Then, V^(x) = V*(x) for all x. 



To end the section, we give the proofs of Theorems 4.1 and 4.2. 

Proof of Theorem 4.1. The condition linx^oo h'(x) = oo implies that b* < oo. Clearly, 
it follows from the definition of V&* and Lemma 4.2 that for x > 

{T-a)V b *{x) < 0, 
(l-^(x))(r-a)V6.(a:) = 0, 
VUx) > 1, 

which in turn imply that for all x > 

max{(r - a)V b *{x), 1 - V^{x)} = 0. 
Hence, the result is a direct consequence of Lemma 4.3. □ 

Proof of Theorem 4.2. If 7r is completely monotone, then ip is also completely mono- 
tone. Hence, ip(x) admits the following representation 

POO 

4>(x) = / e- sx n(ds), (4.3) 
Jo 

where [i is a Borel measure on [0, oo). Making use of Cramer's estimate for ruin probability 
(see Lemma 2.2), we have 



${x) ~ Ce- {R+p(5))x , x 



oo, 



where R is a positive constant such that \l/(— R) = 0, and C is a nonnegative constant. 
Consequently, we have 

lim e p{5)x i){x) = 0. 



x— >co 

P(S) 



This together with (14.31) give J Q P n(ds) = 0, and hence 



oo 



e PW x $(x) = / e~ sx+p{&)x ii(ds) = / e~ tx fi{p{5) + dt), 

Jp(8) JO 

which is completely monotone. In particular, we have (e p ^ x ip(x))'" < 0. Thus, 
h"'(x) = p(5) 3 e p(5)x - {e p{5)x i!{x))"' > 0, x > 0. 

That is, h' is strictly convex on (0,oo). The result follows from Theorem 4.1 since 
lim^oo h'(x) = oo. □ 
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5 Examples 



In this section, we present an example with mixed-exponential jump-diffusion process and 
a list of completely monotone Levy densities. 

Example 5.1. (Mixed-exponential jump- diffusion process - see also Asmussen et al. 
(2004) and Mordecki (2004)) 

Consider the process X = {X t : t > 0} given by 

X t = at + aB t + Y^Yki (5.1) 

k=l 

where B = {B t : t > 0} is a standard Brownian motion, N = {N t : t > 0} is a Poisson 
process with parameter A, Y — {Y k : k > 1} is a sequence of independent and identically 
distributed random variables with density 

<x) J>-zu*mt*- *>«. (5 . 2) 

with p, q > 0, p + q — 1, r/^, ^ > 0, Aj, Bj > 0, J2jLi Aj = 1, 5^" =1 -Bj = 1- As usual, 
we assume that the processes B, N, and Y are independent. When m = n = 1, the 
process A reduces to the double-exponential jump- diffusion process (see Kou and Wang 
(2003)). Since 7r defined in (15. 2p is completely monotone, a direct application of Theorem 
4.2 shows that the barrier strategy at b* is an optimal strategy. □ 

Besides Example 5.1, other examples of Levy processes with completely monotone 
densities can be found in the literature. Several of them are listed below (for details, 
see Bagnoli and Bergstrom (2005), Loeffen (2008), Yin and Wang (2009) and Kyprianou, 
Rivero and Song (2010)): 



• a-stable process with Levy density: tt(x) = Xx 1 a , x > with A > and a G 
(0,1)U(1,2); 

• one-sided tempered stable process (particular cases include gamma process (a = 0) 
and inverse Gaussian process (a = 1/2)) with Levy density: tt(x) = \x~ 1 ~ a e~ /3x , x > 
with p, A > and -1 < a < 2; 

• the associated parent process with Levy density: tt(x) = \ix~ 1 ~ a e~ l3x +\2X~~ 2 ~ a e~ f5x , 
x > with Ai, A 2 > and — 1 < a < 1. 

Note that they all satisfy the condition J °° 7i(x)dx = oo. Moreover, some distributions 
with completely monotone density functions are given below: 

• Weibull distribution with density: f(x) = crx r ~ 1 e~ cxr , x > 0, with c > and 

< r < 1; 
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Pareto distribution with density: f(x) = a(l + x) a 1 , x > 0, with a > 0; 

mixture of exponential densities: f(x) = ^2™ =1 Aif3ie~ l3iX , x > 0, with Ai > 0, > 
for i = 1, 2 • • • , n, and X^=i ^ = ^ 



• gamma distribution with density: 



x c-l e -x/P 



fix) = — . , X > 0, 

JK ' r( c )/3 c ' 



with /3 > and < c < 1. 



Remark. For a spectrally negative Levy process, whenever the Levy measure has a 
density which is log-convex, then Kyprianou et al. (2010) showed that the optimal strategy 
is a barrier strategy. Since there is no condition on the upward jumps in Theorem 4.2, 
it might be possible to generalize the result of Kyprianou et al. (2010) to the present 
situation. That is, for any Levy process with arbitrary positive jumps, if the Levy density 
of negative jumps is log-convex, then the optimal dividend strategy is a barrier strategy. 
However, we are not able to give a formal proof of such a conjecture. □ 
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